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Abstract. Occurrences of jet engine power loss and dam- 
age have been associated with flight through fully glaciated 
deep convection at —10 to — 50 °C. Power loss events com- 
monly occur during flight through radar reflectivity (Z e ) less 
than 20-30 dBZ and no more than moderate turbulence, of- 
ten overlying moderate to heavy rain near the surface. During 
2010-2012, Airbus carried out flight tests seeking to charac- 
terize the highest ice water content (IWC) in such low-Z e 
regions of large, cold-topped storm systems in the vicinity 
of Cayenne, Darwin, and Santiago. Within the highest IWC 
regions encountered, at typical sampling elevations (circa 
1 1 km), the measured ice size distributions exhibit a notably 
narrow concentration of mass over area-equivalent diame- 
ters of 100-500 pm. Given substantial and poorly quantified 
measurement uncertainties, here we evaluate the consistency 
of the Airbus in situ measurements with ground-based pro- 
filing radar observations obtained under quasi-steady, heavy 
stratiform rain conditions in one of the Airbus-sampled loca- 
tions. We find that profiler-observed radar reflectivities and 
mean Doppler velocities at Airbus sampling temperatures are 
generally consistent with those calculated from in situ size- 
distribution measurements. We also find that column simula- 
tions using the in situ size distributions as an upper boundary 
condition are generally consistent with observed profiles of 
Z e , mean Doppler velocity (MDV), and retrieved rain rate. 
The results of these consistency checks motivate an examina- 


tion of the microphysical pathways that could be responsible 
for the observed size-distribution features in Ackerman et al. 
(2015). 


1 Introduction 

Between 1990 and 1998, commuter transport jets experi- 
enced at least 10 incidents of engine power loss near thunder- 

storms at elevations of 8.5-9.5kma.m.s.l. and temperatures 
of —27 to —33 °C (Lawson et al., 1998; Mason et al., 2006). 
During these incidents the engines lost power slowly prior 
to abrupt uncommanded reduction of power to idle, com- 
monly referred to as rollback, and in most cases engine power 

authority was automatically returned after descent through 
the melting level. Crews reported no more than trace oc- 
currences of hail or graupel. In 1997, industry-sponsored 

flight tests successfully reproduced a rollback event while 
circling downwind of continental deep convection at 8.8 km 
and —28 °C over Alabama, induced by enhanced ice wa- 
ter content (IWC), sustained at likely 0.5-1 gm~ 3 or more 
over a 10 km length scale, based on estimates obtained in- 
directly owing to primary probe failure (Strapp et al., 1999; 
Mason et al., 2006). On the one hand, the flight tests were 
deemed a success since a second, modified engine simulta- 
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neously experienced no rollback, and no further events have 
been reported on that particular aircraft since such modifica- 
tions were made standard (Mason et al., 2006). On the other 
hand, the industry concluded (e.g., Dezitter et al., 2013) that 
the inability of current instruments to reliably and accurately 
measure such high IWC posed a major obstacle to addressing 
similar hazards being identified in other aircraft types (Strapp 
et al., 1999, 2005; Mason et al., 2006). 

After a period of confusion regarding the potential role of 
supercooled liquid water owing to the appearance of stream- 
ing meltwater on heated cockpit windshields during some 
events, by 2003 the connection between both commuter and 
large transport jet engine power loss events and the inges- 
tion of copious ice under glaciated conditions was estab- 
lished, ultimately involving a dozen engine types (Mason and 
Grzych, 201 1). As of 2011, a database of roughly 100 docu- 
mented events of jet engine power loss or damage connected 
with the engine core had indicated several patterns in the 
meteorological conditions encountered (Grzych and Mason, 
2010; Mason and Grzych, 201 1). Namely, events occurred al- 
most exclusively near deep convection, under hypothesized 
and sometimes unambiguously fully glaciated conditions, 
in the absence of pilot-reported equivalent radar reflectiv- 
ity (Z e , less than 20-30 dBZ) or more than moderate turbu- 
lence at flight level. Hereafter we will refer to event-related 
microphysical conditions as “high IWC-low-Z e ” conditions, 
where the IWC designation remains qualitative or relative, 
as discussed further below, and the Z e threshold is based on 
radar sensitivity, without any specified relation to convective 
or stratiform structures 

Outside of these commonalities, events occupied a rather 
wide envelope of storm conditions: 38% occurred within 
an oceanic mesoscale convective system (MCS), where pi- 
lots typically flew directly through the central and deepest 
clouds owing to a complete lack of detectable Z e at flight 
altitude and no identification otherwise of flight-level con- 
vective cores; 34% occurred within a strong tropical MCS, 
where pilots detected convective cores at flight altitude, but 
maneuvered to avoid them, typically passing over regions of 
heavy rain; 8 % occurred in a continental MCS, where pilots 
again detected and diverted around cores at flight altitude, 
flying over heavy rain regions; 6 % occurred in classical con- 
tinental anvils, where cores were identified and avoided on 
the downwind side of storms, in this case without significant 
rain immediately below; and the remainder occurred vari- 
ously within tropical multicell, elevated warm frontal, vigor- 
ous winter lake-effect, or developing, immature, continental 
convection (Mason and Grzych, 201 1). 

To evaluate event distribution geographically and season- 
ally, Grzych and Mason (2010) derived maps of high IWC- 
low-Ze frequency, based on statistical comparison of reanal- 
ysis model fields with the following event conditions: high 
precipitable water (median of 6 cm), marginal to modest at- 
mospheric instability (median convectively available poten- 
tial energy of 1400Jkg _1 ), and modest wind shear (gener- 


ally sufficient to support squall lines). When also consider- 
ing flight traffic frequency, these maps generally explained 
the concentration of events within 45° of the equator under 
preferentially moist conditions seasonally. 

It is notable that the MCS class events documented (80 % 
of the total) occurred beneath a cirrus anvil shield of at least 
185 km in maximum dimension. However, an analysis of 
temperature instrument errors caused by high IWC-low-Z e 
conditions during some events has suggested that enhanced 
IWC regions might commonly be as short as (9(1 0 km) in 
length, making it impossible to rule out the sufficiency of 
either long or short duration exposure to cause such events 
(Mason et al., 2006). An exhaustive analysis of conditions 
encountered with and without power loss events led to only 
one consistently recommended hazard detection action avail- 
able using current cockpit instrumentation: use the aircraft 
radar’s tilt and gain to scan below the freezing level for mod- 
erate to heavy rain (Grzych and Mason, 2010; Mason and 
Grzych, 2011). 

2 Microphysical conditions 

From the perspective of industry, the most recent effort to 
rigorously establish the condensed water content possible in 
convection was undertaken by the Royal Aircraft Establish- 
ment in the 1950s using a pitot-type tube that collected ice 
and water condensate in a thermostatically controlled heated 
tube and delivered it to a measuring system inside the air- 
craft (McNaughton, 1959; Mason et al., 2006). Measure- 
ments were made at roughly 10 km resolution over a tem- 
perature range of 0 to —26 °C, with little reference to Z e or 
degree of glaciation encountered. Regions of the highest ex- 
pected turbulence and water content were avoided for rea- 
sons of safety, but the flight paths were also adjusted to find 
and maintain high water content in order to challenge engine 
performance. Based on roughly 200 cloud traverses in the 
vicinity of Entebbe, Singapore, and Darwin, maximum total 
water contents encountered were in excess of 6 and 5 gm -3 
over length scales of 10 and 100km, respectively, and more 
than 5 % of all reported measurements at each location ex- 
ceeded 2gm -3 , with an accuracy of 15% estimated at that 
time. 

Motivated by the 1990s commuter jet engine power 
loss problem, Lawson et al. (1998) reanalyzed more lim- 
ited in situ measurements collected within anvils over 
Montana during the Cooperative-Convective-Precipitation- 
Experiment (CCOPE) (Knight, 1982) and near the Tiwi 
Islands during the Central Equatorial Pacific Experiment 
(CEPEX) (Heymsfield and McFarquhar, 1996; McFarquhar 
and Heymsfield, 1996). Since no direct measurements of 
IWC were obtained during these campaigns, it was calcu- 
lated from the most suitable data available. From CCOPE, 
IWC was calculated from a particle measuring systems 
(PMS) two-dimensional precipitation (2DP) probe and scan- 
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ning radar reflectivity fields co-located with flight tracks, 
with an estimated accuracy of 50 % (Heymsfield and Palmer, 
1986; Heymsfield, 1986). From CEPEX, IWC was calculated 
from a PMS two-dimensional cloud particle probe (2DC), 
with an estimated accuracy of a factor of 2, owing largely 
to sensitivity to mass-dimension relationships applied, plus 
roughly 50 % attributed to particle counting uncertainty (Mc- 
Farquhar and Heymsfield, 1996). The 2DP and 2DC probes 
respectively measure particles that are nominally 300-9600 
and 30-1 800 pm in maximum dimension of randomly ori- 
ented projected area (D max ; the diameter of the smallest cir- 
cle enclosing two-dimensional projections of ice particles 
imaged at random orientation in situ). Working in the ab- 
sence of recorded Z e at flight level from these scientific re- 
search campaigns, Lawson et al. (1998) reported peak anvil 
IWC of 1-3 gm -3 , with a suggested factor-of-2 uncertainty, 
falling to well below 1 gm' 3 at a distance of 50 km from 
storm locations of highest Z e , a vicinity commonly reported 
for initiation of commuter and large transport jet engine 
events (Mason et al., 2006). Ice size -distribution properties 
were not extensively analyzed, except to report the largest 
particles encountered, which were commonly on the order of 
2 mm in D max . 

In a more detailed analysis of ice properties in three 
CEPEX anvils at elevations of 7-14 km and temperatures of 
— 15 to —65 °C, where convection was approached as closely 
as the pilots considered safe, McFarquhar and Heymsfield 
(1996) reported a predominance of irregular, plate-like, com- 
pact spatial crystals and irregular aggregates in the highest 
IWC regions (0.1-1 gm“ 3 , calculated as described above). 
They reported typical mass median D max (MMD max ; the 
D max above and below which half of mass resides) de- 
creasing from 500 to 100 pm with increasing height. Promi- 
nent peaks in the mass size distribution were found at 150- 
600pm. Even if IWC were as high as 4gm -3 , Fig. 1 il- 
lustrates that such relatively small mass median diameters 
would be associated with Z e less than 20dBZ if present 
in monomodal size distributions with typical assumptions 
regarding ice crystal properties (furthermore, we note that 
such calculations do depend strongly on assumed ice crystal 
mass-dimension characteristics, which are not directly mea- 
sured by any current instruments, as well as size distribution 
properties). 

Based on measurements during the Tropical Composition, 
Cloud and Climate Coupling Experiment (TC4) campaign 
over Central America (Toon et al., 2010), reportedly the 
most extensive tropical anvil measurement campaign since 
CEPEX (Lawson et al., 2010), and the NASA African Mon- 
soon Multidisciplinary Analyses (NAMMA) campaign over 
West Africa (http://airbornescience.nsstc.nasa.gov/namma/), 
the ice mass size distributions in fresh anvils, aged anvils, 
and convective turrets alike were consistently found to be 
dominated by crystals with /) lnax of 100-400 pm (Lawson 
et al., 2010). During both TC4 and NAMMA, fresh anvils 
contained columns, capped columns, relatively compact ir- 



Figure 1. For an ice water content of 4 g m~ 3 with a monomodal ice 
number size distribution that is either monodisperse or exponential 
in shape, equivalent radar reflectivity (Z e ) as a function of mass 
median maximum dimension (MMD max ) for ice particle properties 
following Baker and Lawson (2006, their Table 1) and mass me- 
dian area-equivalent diameter (MMD e q) for ice particle properties 
following Eq. (1) (Locatelli and Hobbs, 1974). 


regular crystals and irregular aggregates that were gener- 
ally smaller than 2 mm in Z) max , whereas turrets contained 
1-2 mm graupel over a wide temperature range (—10 to 
—50 °C). Considering this diversity of new data with strik- 
ingly similar microphysical features, Lawson et al. (2010) 
attributed the predominance of several hundred micron par- 
ticles to cloud droplets freezing heterogeneously at temper- 
atures colder than — 12 °C and then growing by vapor de- 
position, while smaller and larger outflow ice preferentially 
evaporated and sedimented, respectively. 

During TC4, IWC was calculated by applying the habit- 
independent area-dimensional relation derived in Baker 
and Lawson (2006) to measurements from an underwing 
mounted PMS two-dimensional stereo (2DS) probe over 
a size range of roughly 10-3000 pm in Z) max . Calculated 
IWC was shown to agree with direct measurements from 
a counterflow virtual impactor (CVI) inlet mounted on the 
fuselage of the aircraft to within 20% on average at IWC 
concentrations below 0.5 gm -3 , the estimated saturation of 
that instrument under TC4 measurement conditions (Law- 
son et al., 2010). Within a heated CVI inlet, impacted ice 
is evaporated into a dry airstream, which is in turn sampled 
by a Lyman-a hygrometer; IWC uncertainty is estimated to 
be 10 % for values larger than 0.2 gm -3 (Twohy et al., 1997). 
Lawson et al. (2010) do not offer an estimate of uncertainty 
in IWC calculated or measured during TC4 or NAMMA. Al- 
though IWC exceeding 2 g rrT 3 was calculated in one con- 
vective turret, maxima were less than 1 gm -3 in fresh and 
aged anvil regions. 

In a data set compiled for development of satellite retrieval 
algorithms, comprising thousands of 5 s averaged aircraft 
data from TC4, NAMMA, and other campaigns, IWC values 
greater than 2 gm -3 are very rare (Delanoe et al., 2014). Fur- 
thermore, we note that no research aircraft involved encoun- 
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tered engine power loss events although high IWC-low-Z e 
conditions were not specifically avoided, to our knowledge. 

Comparison of IWC retrievals using CloudSat data 
(25° S-25° N) from 7 July to 16 August 2006 over the tropics 
(25° S-25°N; ~ 1.8 km x 1.4 km x 0.5 km resolution) showed 
rare occurrences of greater than 2.5 gm~ 3 at 10km in eleva- 
tion (frequency of less than 10 -8 ) from two out of three algo- 
rithms (Wu et al., 2009). However, a third algorithm showed 
occurrences of IWC exceeding 10 gm , and the spread of 
algorithm results diverged sharply at IWC values greater than 
1 gm -3 . 

Motivated directly by the engine power loss problem, 
a more recent satellite-based study aimed to establish IWC 
conditions conceivably encountered by commercial aircraft 
flying at an elevation of 10 km in the tropics using Cloud- 
Sat and CALIPSO measurements and retrievals (Gayet et al., 
2014). After identifying and analyzing a well-observed trop- 
ical MCS off the equatorial coast of Brazil, they concluded it 
likely that the storm included an area wider than 55 km with 
IWC greater than 1 gm -3 and 94GHz Z e below 18dBZ, in- 
cluding 2.5 km footprint peaks of 24gm -3 , with a large un- 
certainty owing to unknown ice morphological and size dis- 
tribution properties. For an ice mass mixing ratio of 4 gm“ 3 . 
Fig. 1 illustrates the relationship of MMD max to reflectiv- 
ity for the simplest case of monodisperse ice particles with 
mass following Eq. (1), and for exponential size distributions 
with slopes over an observation-based range of 1-1000 cm -1 
(Heymsfield et al., 2013). Holding size distribution proper- 
ties fixed, denser ice particles would lead to increasing re- 
flectivity for the same MMD max and ice mass mixing ratio. 

Focusing briefly on observed ice morphology, the McFar- 
quhar and Heymsfield (1996) and Lawson et al. (2010) stud- 
ies both identify a sharp distinction between crystal shapes 
seen in deep convection outflow vs. those seen in cirrus that 
are formed in situ. Although both may contain a mixture of 
habits that include irregular and rather compact appearing 
particles of several hundred microns in D max , cirrus formed 
in situ are most commonly characterized by a notably larger 
fraction of rosettes (Lawson et al., 2010, per their Fig. 2). Al- 
though not reported from TC4 or NAMMA, deep convection 
that is electrically active may also be distinguished by a pre- 
dominance of chain aggregates (Stith et al., 2004), which 
have been found to be composed of small frozen droplets or 
faceted crystals with D max commonly less than ~ 25-100 pm 
(e.g., Connolly et al., 2005; Gayet et al., 2014; Stith et al., 
2014). The occurrence of chain aggregates is attributed to 
the presence of electric fields sufficient to enhance the ag- 
gregation process, as well as sufficiently many such small 
crystals to aggregate, perhaps owing to homogeneous freez- 
ing of abundant cloud droplets in the strong, mixed-phase up- 
drafts associated with electrical activity (e.g., Connolly et al., 
2005). As discussed further below, the complexity of crystal 
morphology contributes substantially to uncertainty in the re- 
lationship between IWC, ice mass size distribution, and Z e in 
clouds associated with deep convection. 


In this two-part study, we examine a set of measure- 
ments gathered by Airbus during 2010-2012 in deep convec- 
tion near Cayenne, Darwin, and Santiago. As described by 
Grandin et al. (2014), an Airbus objective was to obtain pre- 
liminary measurements to evaluate newly proposed standards 
for engine performance in glaciated icing conditions (Maz- 
zawy et al., 2007). Conditions sampled included oceanic and 
continental conditions, monsoonal and sea breeze convec- 
tion, and a range of storm size extents. Similar to the Mc- 
Naughton (1959) campaigns, the purpose was to find and 
stay within the highest IWC regions possible. Here in Part 1, 
we first survey the IWC, mass size distribution, and Z e de- 
rived from measurements; these are the factors immediately 
relevant to characterizing high IWC-low-Z e conditions. We 
compare the flight test data with previous literature and with 
remote-sensing measurements of a well-observed MCS over 
Darwin during the Tropical Warm Pool International Cloud 
Experiment (TWP-ICE) campaign (May et al., 2008). Fi- 
nally, we use Airbus observations to initialize steady-state 
column simulations of heavy stratiform rain columns, and 
compare predicted Z e and Doppler fall speed profiles with 
observations. Aside, we note that although many power loss 
events may be within a region classified as stratiform rain 
by some measure (Grzych and Mason, 2010), it is not clear 
whether power loss events are associated with classical strat- 
iform rain columns (Biggerstaff and Houze, 1991), com- 
monly associated with a radar reflectivity bright band (e.g., 
Steiner et al., 1995), or some type of transitional regime. 
A transitional regime region is already difficult to objec- 
tively define in a mid-latitude squall line (cf. Biggerstaff and 
Houze, 1993), where convective and stratiform rain regions 
are commonly well separated, and is more complicated in 
large tropical systems owing to the variety of convective re- 
gion structure (e.g., Mapes and Houze, 1992). Here simu- 
lations are limited to heavy stratiform rain regions in order 
to make use of Doppler fall speed information under rele- 
vant quiescent conditions where mean vertical wind can be 
neglected at elevations circa 1 1 km relative to reflectivity- 
weighted ice fall speed. 

3 Airbus in situ measurements 

The Airbus flight campaigns included two in situ probes in- 
tended to maintain functionality under high IWC conditions: 
a newly designed Science Engineering Associates (SEA) Ro- 
bust solid hot-wire probe used for estimation of ice mass 
mixing ratio (cf. Strapp et al., 2005) and an Airbus imag- 
ing nephelometer designed for measurement of liquid droplet 
size distribution but used here for estimation of ice mass size 
distribution (Roques, 2007; Dezitter et al., 2013; Grandin 
et al., 2014). 

Since about 2005, it has been realized that hot-wire de- 
vices using trapping geometries are subject to relatively 
low capture efficiencies owing to re-entrainment of only 
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a fraction of particles that bounce out of the capture vol- 
ume (Emery et al., 2004; Strapp et al., 2005; Korolev et al., 
2013b). The Robust probe was developed primarily to pro- 
vide surrogate IWC measurements for the calibration of 
a wind tunnel producing high IWC at relatively high wind 
speed (Strapp et al., 2008), conditions where essentially all 
available hot-wire probes were known or suspected to fail 
owing to particle impact damage (Strapp et al., 1999). In 
a wind tunnel. Robust probe measurements were intended 
to be used alongside ice ingestion measurements to deduce 
absolute IWC and ice mass size distributions. During the 
course of calibrations, the efficiency of the capture, melting, 
and evaporation by the Robust probe for shaved ice parti- 
cles was estimated at about 0.4, the value assumed here. The 
probe was subsequently recommended to Airbus on the basis 
of durability and its simple principle of operation, requiring 
a relatively simple installation, to provide real-time guidance 
during flight tests. 

The accuracy of applying a 0.4 collection efficiency to es- 
timate IWC from the Robust probe under natural ice particle 
conditions remains unknown and will require further inves- 
tigation (Grandin et al., 2014). Comparisons with IWC esti- 
mated from the imaging nephelometer in Airbus flight tests 
has indicated smaller efficiencies of 0. 2-0.3 (Grandin et al., 
2014). Although the IWC accuracy was initially estimated as 
30 % (assuming the primary source of error to be the loss of 
re-entrained mass, such that measurements uncorrected for 
efficiency should provide an accurate lower limit on IWC), 
we can offer no further evidence of accuracy here to establish 
anything better than the raw factor of 2 uncertainty owing to 
overall collection efficiency. Recent flight tests have provided 
a data set that is hoped will improve efficiency characteriza- 
tion in the near future (Grandin et al., 2014). Because the 
probe’s efficiency for liquid water content is close to unity, 
such that liquid and ice contributions to mass cannot be ac- 
curately determined, in Fig. 2 we include only measurements 
taken at temperatures colder than — 20 °C, where airframe 
icing was non-existent or negligible; size-distribution mea- 
surements are limited to temperatures colder than — 33°C, as 
discussed further below. 

In Fig. 2 the total condensed water content (CWC; liq- 
uid plus ice) reported by the Robust probe during each cam- 
paign is plotted as a function of elevation over three averag- 
ing lengths: 0.2 km (1 Hz data as reported, at typical airspeed 
of 200 ms -1 ), 2.5, and 10km. Data include eight flights from 
Cayenne (F1419-F1426), six flights from Darwin (F1403- 
F1409), and nine flights from Santiago (F1539-F1547). We 
note that the Cayenne and Darwin flights targeted cruise el- 
evations circa 10 km, whereas the Santiago flights targeted 
warmer temperatures; the concentration of IWC as a function 
of elevation at each location should therefore be viewed pri- 
marily as an indication of sample frequency. Based on a more 
detailed analysis of the Cayenne and Darwin Robust probe 
data as a continuous function of length scale within iden- 
tified legs, Grandin et al. (2014) report that 99th percentile 


Cayenne 0.2-km Darwin 0.2-km Santiago 0.2-km 



02468 02468 02468 

CWC (g nr?) CWC (g nr) CWC (g nr) 

Cayenne 2.5-km Darwin 2.5-km Santiago 2.5-km 



02468 02468 02468 

CWC (g rn ; CWC (g nr) CWC (g nr s ) 
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CWC (g m ’) CWC (g m ) CWC (g in' 3 ) 


Figure 2. Total condensed water content (CWC) measurements by 
the Robust probe during flight tests from Cayenne (left column), 
Darwin (middle column), and Santiago (right column). Measure- 
ments are included only for air temperatures were below — 20 °C, 
where liquid water contributions are expected to be negligible. Con- 
sidering a typical air speed of 200 ms -1 , 1 Hz measurements cor- 
responding to a horizontal length scale of roughly 0.2 km (top row) 
slightly exceed those averaged over 2.5 and 10 km length scales. 


CWC is a surprisingly weak function of length scale over 2 
orders of magnitude (0.2-20 km), consistent with the pattern 
seen here across all three locations in Fig. 2. They also show 
that a leg with statistically extreme peak CWC does exhibit 
a sharper decrease in CWC with averaging distance, as might 
be expected within a narrow updraft, for instance. 

Taken together, the envelope of measurements is quali- 
tatively similar to that reported by McNaughton (1959) in- 
sofar as 10 km mean condensate mixing ratios in excess of 
2gm -3 are found in all three locations. At a length scale of 
10km, 2gm -3 is roughly twice the maximum value shown 
across the three CEPEX anvils examined by McFarquhar and 
Heymsfield (1996, cf. their Fig. 3). During TC4, 2 gm -3 was 
only reported at a length scale of 0(1 km) in an updraft ex- 
ceeding 20ms -1 (Fawson et al., 2010). TC4 IWC at a 10km 
length scale reached at least 0.25 gm -3 in fresh anvils, and 
mean anvil IWCs reported appear generally consistent with 
CEPEX data (McFarquhar and Heymsfield, 1996; Fawson 
et al., 2010). On the face of it, these results may seem statis- 
tically inconsistent; however, at least some degree of differ- 
ence could be explained by the objective of the Airbus flight 
tests to sample very large, cold-topped storms and to pri- 
oritize the repeated measurement of the highest IWC found 
there. A large database of such measurements does not cur- 
rently exist. Other differences may be within the large and 
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Figure 3. Cayenne flight F1423 track over brightness temperature from Meteosat 9 channel 10 at 12:45 UTC on 16 May 2010. 


poorly quantified uncertainty of both measured and calcu- 
lated IWC values. 

We next consider in greater detail the measurements ob- 
tained during a flight that encountered relatively high IWC. 
From the track of Cayenne flight F1423 (Fig. 3), with spi- 
rals in an identified high IWC-low-Z e region at circa 11 km 
in altitude. Fig. 4 shows ice particle number size distribu- 
tions derived from the nephelometer measurements. At tem- 
peratures circa —43 °C during the flight segment analyzed 
(Fig. 5), the sampled cloud is fully glaciated. While capped 
columns are commonly present (see Ackerman et al., 2015, 
hereafter Part 2), the majority of crystals appear irregular, 
as found in CEPEX anvils (cf. McFarquhar and Heymsfield, 
1996), and the nephelometer images do not commonly pro- 
duce images of sufficient clarity to distinguish rime or other 
morphological details. The nephelometer size distributions 
are reported as a function of area-equivalent diameter (D e q ), 
defined as the diameter of a circle with area equal to the two- 
dimensional projection of an ice particle imaged at random 
orientation in situ, over a bin midpoint range of 67.5 pm to 
3 mm. Also shown is the mass size distribution calculated as- 
suming a widely used relationship of particle m in mg to D eq 
in mm derived from measurements of natural crystals sam- 
pled in the Cascade Mountains, Washington (Locatelli and 
Hobbs, 1974): 

m = 0.037 D^ 9 . (1) 


Here we apply this relation using the same dimension 
D eq defined by Locatelli and Hobbs (1974) for their single- 
particle measurements, applicable both to aggregates of 
densely rimed radiating assembles of dendrites or dendrites 
(£> e q of 2-12 mm measured) and to aggregates of unrimed 
radiating assemblages of plates, side planes, bullets, and 
columns ( D eq of 1-3 mm measured). Aside, we note that this 
D eq is slightly different from the average of chord lengths 
parallel and perpendicular to a photodiode array substituted 
in Brown and Francis (1995). Whereas this m-D eq relation 
was derived from particles with D eq of 1-12 mm (based on 
roughly fifty particles total), it is applied here at all D eq , with 
the exception that effective density is limited to that of bulk 
ice, which is relevant where D eq is smaller than ~ 100 pm; 
such small particles do not contribute substantially to the 
mass size distributions shown here, though bin-wise uncer- 
tainty in mass size distributions calculated using this ap- 
proach is not quantifiable at this time. As discussed above, 
roughly a factor of 2 uncertainty in calculated IWC may be 
associated with uncertainty in the validity of Eq. (1) or an- 
other such relationship (e.g., McFarquhar and Heymsfield, 
1996). Also shown in Fig. 4 are Rayleigh radar reflectivities 
calculated from the ice mass size distributions as the inte- 
gral of the sixth moment of the melted diameter multiplied 
by a dielectric factor of 0.208/0.93 (Smith, 1984). 

As mentioned above, here we focus primarily on the mass 
and radar reflectivity that contribute to high IWC-low-Z e 
conditions. Figure 5 shows the time series of integrated 
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Figure 4. Cayenne flight F1423 ice number size distribution (top 
row), mass size distribution (middle row), and equivalent radar re- 
flectivity size distribution (bottom row) in three ice water content 
(IWC) ranges: 1-2 gm -3 (left column), 2-4 gm -3 (middle col- 
umn), and > 4gm -3 (right column). Measurements from the an- 
alyzed time period 11:30-13:20 UTC on 16 May 2010, during spi- 
rals within an MCS (cf. Figs. 3 and 5). Red line indicates bin-wise 
median size distribution. 


IWC, mass median Z5 eq (MMD eq ), and Z e obtained from the 
measured size distributions. Here it is shown that IWC ob- 
tained from the Robust probe closely tracks IWC obtained 
from integration of calculated mass size distributions over 
roughly 0. 1 —4 g rrT 3 reported, with the mass calculated from 
the nephelometer exceeding that from the Robust probe by 
roughly 25 % at the highest values reported. Calculated Z e 
values are below 15 dBZ, although the sharp drop in Z e size 
distributions at the largest particles sizes indicates that the 
nephelometer size distributions do not include the largest par- 
ticles contributing to calculated Z e in all three mass ranges 
shown, so this calculated Z e is low by an unknown amount. 
However, mass size distributions do not appear to be as trun- 
cated as Z e size distributions; although they must also be 
truncated to some degree by an unknown amount, agreement 
of integrated IWC with the Robust probe suggests that this 
may not be a great amount. We might expect a greater mass 
contribution from unsampled larger particles at lower eleva- 
tions than sampled during this flight. 

The most striking feature of the median size distribu- 
tions shown is the consistent peak in particle mass at Z) eq 
of roughly 250 pm in Fig. 4. Although the peak is narrower 
at higher IWC, the time series in Fig. 5 show that calcu- 
lated MMD eq scarcely deviates from 200 to 300 pm over 






11.5 12.0 





Figure 5. Static air temperature, ice water content (IWC), equiva- 
lent radar reflectivity (Z e ), ice mass median area-equivalent diam- 
eter (MMD eq ), and Z e -weighted fall speed (Vz) during Cayenne 
flight F1423. 


more than an order of magnitude variability in IWC. Fig- 
ure 6 further demonstrates that when the calculated mass 
size distributions are averaged within the three mass ranges 
shown in Fig. 4 (1-2, 2 — 4, and greater than 4gm -3 ) and 
then normalized by their total mass for Flight 1423, they ap- 
pear roughly self-similar insofar as they exhibit similar width 
and shape across roughly a factor of 4 in mass concentra- 
tion. Three other flights with the highest IWC encountered 
in each flight test location show similar behavior, with cal- 
culated MMD eq similar to 250 pm, although more smaller 
and larger particles appear likely to have been present in the 
Santiago flight, which sampled lower elevations (see Fig. 2); 
owing also to the likelihood that unsampled larger particles 
contributed more substantially to that flight, Santiago results 
are omitted from Part 2 of this study. Although Robust probe 
data are available at warmer temperatures (see Fig. 2), we 
note that Cayenne and Darwin size distributions are avail- 
able only for flight legs at temperatures colder than — 40°C 
(e.g.. Fig. 5), aside from several short segments of flight 1422 
reaching —33 °C, which are not qualitatively different. 

Within any single flight through dense anvil at circa 1 1 km 
altitude in the Cayenne and Darwin flights, a stable self- 
similarity of mass size distributions over a wide range of 
IWC can be considered a surprising result because it is gen- 
erally thought that the process of ice outflow from convective 
updrafts and subsequent sedimentation consistently leads to 
preferentially larger particles and more mass occurring closer 
to updraft sources (e.g., Lawson et al., 1998, 2010). Consis- 
tent with that reasoning, McFarquhar and Heymsfield (1996) 
reported that larger ice crystals were found preferentially 
closer to convective cores and lower in the anvil, consis- 
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Figure 6. Ice mass size distributions in three ice water content 
(IWC) ranges (1-2, 2^1 and > 4gm -3 shown in cyan, blue and 
red, respectively) based on bin-wise median measurements during 
four flights. Mass concentrations in gm” 3 (left panels) and normal- 
ized by IWC (right panels). Mass median area-equivalent diameter 
(MMD e q) listed for the highest concentration range (> 4gm -3 ). 


Encountering a consistency of MMD eq in the densest anvil 
clouds in different parts of the world may also be surpris- 
ing given the relatively wide range of updraft strengths ex- 
pected (e.g., Varble et al., 2015). One possible explanation is 
that the Airbus flight tests systematically targeted central re- 
gions of the largest and coldest-topped clouds in all locations 
(in some cases very shortly after electrical activity ceased), 
and that the coupled dynamical and microphysical processes 
that lead to the highest IWC in such clouds do so consis- 
tently by producing a concentration of mass at MMD eq of 
roughly 250 pm. Lawson et al. (2010) also reported a con- 
sistent concentration of mass in the 100-400 pm D max range 
in both fresh and aged anvils around Central America and 
West Africa; a particularly stable peak of mass in that size 
range was also observed within turrets, leading them to sug- 
gest a microphysical pathway of ice formation in updrafts 
that could explain such a pattern. The size range of mass con- 
centration found during the Airbus flight tests is also roughly 
consistent with that reported from CEPEX, although IWC 
was substantially less throughout previous scientific cam- 
paigns, even within graupel-containing updraft turrets sam- 
pled. Leroy et al. (2015) more recently report MMD eq typ- 
ically 250-500 pm and weakly decreasing with increasing 
IWC over 0.5-3 gm -3 at — 36 °C in a system extensively 
sampled near Darwin during the recent High Altitude Ice 
Crystals/High Ice Water Content campaign, but similar mea- 
surements in another system yielding a typical MMD eq of 
400-600 pm instead weakly increasing with IWC. Aside, we 
note that shattering artifacts that may contaminate airborne 
probe measurements are relatively reduced for higher-order 
moments of the size distribution such as mass (Korolev et al., 
2013a; Jackson and McFarquhar, 2014). Since the Korolev 
et al. (2013a) study was performed for probes with different 
inlet configurations, we expect that general conclusion can be 
extended to the Airbus nephelometer. It has been found that 
size-distribution measures such as MMD max may be subject 
to roughly 20 % uncertainty owing to shattering artifacts, for 
instance (Jackson and McFarquhar, 2014). 


4 Profiling radar measurements 


tent with size sorting documented by Garnett et al. (2005) 
within Florida anvils. However, horizontal regions of rel- 
atively uniform MMD max were also noted in all three of 
the CEPEX anvils examined by McFarquhar and Heyms- 
field (1996), leading the authors to conclude that size sorting 
must not be the only dominant process in the CEPEX anvils 
they examined. A later study of precipitating deep convection 
spanning several field campaigns also reported that exponen- 
tial fits to ice size distributions exhibited order-of-magnitude 
variability in intercept but little in slope within individual spi- 
ral loops of 5-10 km diameter (Heymsfield et al., 2002). 


Owing to the uniqueness of the Airbus flight test data, it is de- 
sirable to corroborate it somehow with independent measure- 
ments. It is also desirable to make use of remote-sensing data 
if possible, in part because they commonly offer the advan- 
tage of better spatiotemporal coverage than in situ flight legs, 
and in part because they are free of the various limitations 
associated with in situ probes and measurement techniques. 
Satellite infrared images can offer instantaneous storm wide 
coverage, for instance, but include no direct information 
on either IWC or Z e within dense anvil regions. Ground- 
based scanning weather radars can offer storm wide Z e mea- 
surements at relevant volumetric resolutions of 0(1 km 3 ), 
but variability in ice hydrometeor single-crystal and size- 
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distribution properties preclude reliable retrievals of IWC at 
this time. Vertically pointing radars can also provide non- 
attenuated Z e at relevant volumetric scales of 0 ( 1 00 m 3 ), in 
addition to mean Doppler velocity (MDV). Since MDV rep- 
resents a higher order moment of the ice size distribution than 
Z e , it provides the valuable benefit of an independent mea- 
surement constraint that is also sensitive to mass size distri- 
bution. Here we will focus on extended stratiform regions 
where the mean contribution of vertical winds to MDV can 
be neglected relative to the reflectivity-weighted fall speed of 
ice at elevations circa 1 1 km. The chief disadvantage of pro- 
filing radars is the severely reduced spatial coverage of a pro- 
file compared with a scanning radar volume. However, on 
23 January 2006 during the TWP-ICE campaign, the NOAA 
S-band profiler obtained a time series of measurements dur- 
ing evolution of a tropical MCS near Darwin (Fig. 7). The 
remainder of this work is devoted to examination of these 
measurements in comparison with the Airbus flight test data. 

The S-band MDV curtain in Fig. 8 shows two broad peri- 
ods of stratiform rain passing over the profiler on 23 January 
2006, bisected by a several-hour period of intermittent tur- 
bulence throughout the observed column at roughly 23.65- 
23.83. The stratiform periods are characterized by a sharp 
decrease in MDV at the melting-level height of roughly 5 km, 
where melting ice also creates a bright band in Z e as the 
meltwater undergoes conversion to faster-falling rain. Dur- 
ing the turbulent convective period, Doppler velocities reach 
localized peaks of nearly 10 ms' 1 upward and exceeding 
10ms -1 downward. The strongest drafts are co-located with 
the strongest deep echoes observed in the storm above 15 km. 
Where a bright band is present, owing to the melting of ice 
into stratiform rain, Z e in excess of lOdBZ commonly still 
extends above 1 1 km, with some detectable echoes extending 
above 15 km. Below the bright band, heavy rain is indicated 
by Z e in excess of 30 dBZ. 

Figure 9 shows time series extracted from the S-band pro- 
files and surface rain rate measured by a co-located NOAA 
Joss Waldvogel disdrometer. Reflectivity, MDV, and Doppler 
spectral width are shown at an elevation of 11.7 km, where 
soundings available from Darwin on 23 January indicate 
a temperature of —43 °C, close to that during Cayenne 
flight F1423. Although Z e at this elevation nears 30 dBZ 
(during convective periods), it does not exceed that. The two 
stratiform time periods within dashed lines are the contigu- 
ous periods of at least 15 min where, at 11.7 km, Z e is greater 
than 5 dBZ, MDV is 0-2 ms - 1 , and Doppler spectral width is 
less than 2.5 ms -1 . At 4.7km, Z e is greater than 25 dBZ. The 
beginning of the first period (12:17-14:07 UTC) corresponds 
to the earliest panel in Fig. 7 and the end of the second period 
(20: 1 3-2 1 :49 UTC) to the last panel. Also shown in Fig. 9 are 
retrievals of vertical wind speed obtained from the profilers at 
the highest elevation where retrievals were usually possible 
(9 km) and near-surface rain rate where retrievals can confi- 
dently be made (2.5 km) (Williams and May, 2008; Williams 
and Gage, 2009). Retrieved vertical winds and S-band spec- 



Figure 7. Cloud top temperature derived from MTSAT-1R at 12:00, 
16:33, and 22:33 UTC on 23 January 2006 over Darwin. The darkest 
shade is colder than 200 K, overlying a region that includes the Tiwi 
Islands in each image. 


tral width indicate that the two stratiform periods are quies- 
cent, with minimal turbulence and no strong updrafts. During 
each period, maximum Z e below the melting level exceeds 
40 dBZ and maximum rain rate retrieved at 2.5 km reaches 
nearly 4mmh . At the surface, rain rates measured by dis- 
drometer tend to read sustained peaks circa 2— 4mh -1 near 
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Figure 8. Equivalent radar reflectivity (Z e ) and mean Doppler ve- 
locity curtain measured by the S-band profiler on 23 January 2006. 
Dashed lines indicate duration of two stratiform rain periods defined 
as described in text. 


the midpoints of the stratiform periods rather than at the end- 
points. 

Figure 10 shows MDV vs. Z e in linear coordinates during 
each stratiform period. Mean values of MDV indicate Z e - 
weighted fall speeds of roughly lms^ 1 during each period, 
consistent with those calculated from Cayenne flight F1423 
measurements (Fig. 5) using the method of Heymsfield and 
Westbrook (2010). In the profiler data, coherent oscillations 
in retrieved vertical winds at 9 km and in MDV at 1 1 .7 km are 
likely indicative of gravity waves, and these generate a sig- 
nificant component of scatter around the mean in Fig. 10; ex- 
traction of gravity waves is considered beyond the scope of 
this work. Although it is difficult to draw conclusions owing 
to lack of retrieved winds at 1 1 km, it appears possible that 
MDV may be a weak function of Z e , as suggested by MDV 
calculated from the Airbus data, as expected in the case that 
size distributions are not strongly dependent on IWC. Dur- 
ing the second period in Fig. 10, fall speeds do appear faster 
where reflectivity exceeds 20-30 mm 6 m -3 ; omission of the 
full Z e spectrum apparent in the Airbus measurements, as 
discussed above, may substantially truncate the contribution 
of the largest particles to MDV in the Airbus data. 

Overall, the S-band measurements and retrievals appear 
to generally confirm that regions of low-Z e can be asso- 
ciated with relatively stable MDV of circa lms^ 1 under 
conditions similar to those sampled on flight F1423. MDV, 
which includes air motions, does not vary strongly with 
a factor of 4 change in Z e in these regions, consistent with 
rough self-similarity in size-distribution shape found in the 
in situ Airbus measurements. Figure 9 shows that Z e fre- 
quently exceeds 30 dBZ below these stratiform regions, con- 
sistent with high-Ze conditions beneath low-Z e legs aloft. 
Co-located retrievals of rain rate indicate precipitation ex- 
ceeding 3mmh _1 , consistent with heavy rain expected be- 
neath high IWC conditions. 
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Figure 9. Time series of observations and retrievals from co-located 
S-band, UHF, and VHF profilers on 23 January 2006: radar reflec- 
tivity (Z e ), mean Doppler velocity (MDV; negative downward), and 
Doppler velocity spectral width at 1 1 .7 km; retrieved vertical wind 
speed at 9 km; Z e just below the melting level at 4.7 km; retrieved 
precipitation rate at 2.5 km; and precipitation rate measured at the 
surface. Vertical dashed lines indicate stratiform rain periods as in 
Fig. 8. 
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Figure 10. Mean Doppler velocity vs. equivalent radar reflectivity 
measured by S-band profiler during two stratiform periods shown 
in Figs. 8 and 9. Dashed lines indicate mean values within each 
sample; mean reflectivity reported in dBZ units at lower left. 


5 Column model simulations 

As shown above, the ice size-distribution features found in 
Airbus flight tests are found to be generally consistent with 
other recent in situ measurements in the vicinity of deep 
convection and with profiling radar measurements in a large 
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tropical MCS. Here we take a step further to compare one- 
dimensional simulations of the quasi-steady stratiform rain 
regions observed during TWP-ICE with the profiling radar 
observations. This comparison is intended to allow for an 
evaluation of whether the ice mixing ratios measured in high 
IWC regions are consistent with heavy rain rates retrieved 
below, for instance. As shown in Part 2, three-dimensional 
simulations generally fail to reproduce the observed ice size- 
distribution features in deep convection, with concentration 
of mass narrowly in the 100-500 pm D max range at tempera- 
tures below — 40°C. This exercise therefore also serves to 
test how well a column simulation with size-resolved mi- 
crophysics can perform when initialized with the observed 
size distributions and an observation-based estimate of bin- 
wise ice properties, and when limited to the relatively simple 
dominant processes of sedimentation and aggregation in the 
stratiform region. 

We begin with an atmospheric column with 250 m ver- 
tical grid spacing that is initialized thermodynamically 
with a mean of soundings observed around at Darwin at 
21:00UTC on 23 January 2006, obtained from a variational 
analysis of conditions over the TWP-ICE sounding array 
(Xie et al., 2010). This time is selected to be representative 
of high-Ze, heavy stratiform rain indicated in Figs. 8 and 9. 
Furthermore, we note that no sounding is available at Darwin 
during the second half of 23 January owing to balloon fail- 
ures. The atmospheric column is saturated with respect to ice 
at relevant elevations above 8 km. At the melting-level rela- 
tive humidity is roughly 90 %, increasing above and below 
with respect to ice and water, respectively (not shown). 

We initialize the column with Airbus ice size-distribution 
measurements obtained over Darwin during flight F1405. 
The mass size distributions averaged over regions with re- 
ported IWC greater than 4gm -3 shown in Fig. 6 represent 
a mean IWC of 4.3 gm“ 3 measured between 1 1 and 1 1.5 km. 
The calculated Z e and MDV are shown with blue asterisks in 
Fig. 11. Here we assume that the particles have the mass and 
Z) eq that were shown to agree with the Robust probe mea- 
surements. For simplicity, we assume that the ice particles 
are spherical in shape. Given mass and Z) eq , we therefore 
assume that D max is equal to D eq , the projected area is re- 
lated to Z) eq by definition, and the overall particle aspect ratio 
is unity. Particle fall speeds and pair-wise collision efficien- 
cies are calculated following Bohm (1999), with fall speeds 
modified following Heymsfield and Westbrook (2010) as de- 
scribed by Avramov et al. (201 1). 

For these simulations, we adopt a steady-state simulation 
approach, considering that large and long-lived stratiform re- 
gions are changing slowly over hours in time, therefore ex- 
hibiting roughly equal quantities of mass coming from above 
and exiting below at each level. We neglect large-scale ver- 
tical motion. To speed up simulation time to steady state, 
we initialize all elevations above the melting level with the 
same observed mass size distribution. Vapor exchange with 
hydrometeors is permitted to represent evaporation of ice 




Figure 11. Profiles of reflectivity, mean Doppler velocity (positive 
downward), precipitation rate (liquid equivalent), and condensed 
water content (CWC) observed by the S-band profiler or retrieved 
from S-band, UHF, and VHF radars during the stratiform periods 
identified in Fig. 8 (pink and red) are compared with steady-state 
column model results initialized with size distributions measured 
during Airbus test flights at 11-1 1.5 km over Darwin (blue). Sim- 
ulations use ice-ice sticking efficiencies of 0.01 or 0.005, and the 
baseline sounding or a saturated sounding (see text). 


and liquid in the unsaturated regions, treated as described 
in Fridlind et al. (2012b) and discussed further in Part 2. 
However, we fix temperatures and water vapor mixing ra- 
tios to the observed profile, thereby neglecting cooling at the 
melting level, which would otherwise require application of 
divergence terms to remain stable as observed. Melting of 
ice is treated in a rudimentary fashion, by transferring all ice 
to melted equivalent drops when the column temperature al- 
lows. 

In Fig. 11 simulated Z e , Doppler velocity, precipitation 
rate, and CWC (liquid plus ice) are shown for two values 
of ice-ice sticking efficiency (0.01 and 0.005) and two en- 
vironmental conditions (the observed sounding or a satu- 
rated version) after reaching steady state. Also shown are 
profiles measured by the S-band profiler (Z e and MDV) or 
retrieved with a combination of S-band, UHF, and VHF data 
below the melting level (rain rate and CWC). The observed 
and retrieved data at each elevation represent means of the 
available values over the stratiform time periods shown in 
Figs. 8 and 9. At the top of the simulated columns, Z e and 
MDV represent the steady-state boundary condition. As the 
ice aggregates during sedimentation, the Z e monotonically 
increases. In the simulations, the rate of increase scales with 
the sticking efficiency (£,-), defined here as the probability 
that an ice-ice collision results in aggregation. The probabil- 
ity is taken here to be fixed, although it may be a substantial 
function of environmental conditions and ice morphological 
properties in a manner that is not generally well established 
observationally (e.g., Kienast-Sjogren et al., 2013, and refer- 
ences therein). Owing to the omission of mixed-phase parti- 
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cles from this simplified model, ice instantly melts when the 
temperature is warm enough, and simulations therefore do 
not reproduce the melting band feature in Z e or the smooth 
increase in MDV seen in the observations. 

MDV above the melting level is usually similar to ob- 
served or is underestimated by a few tenths of a meter per 
second. MDV below the melting level is similar to observed 
when Z e is similar to observed, including reproducing a re- 
duction of MDV with distance towards the surface, except 
in one simulation with a saturated sounding, indicating that 
evaporation may substantially increase MDV. However, un- 
derestimation of MDV below the melting level in a simula- 
tion could also be attributable to missing mixed-phase pro- 
cesses (e.g., Phillips et al., 2007), which could enhance ag- 
gregation. 

Simulations are shown using an observed sounding (with 
RH of roughly 90 % at the melting level, as described above) 
and with the same sounding that has been saturated with re- 
spect to ice and water above and below the melting level, re- 
spectively. Simulations with the saturated sounding demon- 
strate that in a steady-state environment without evaporation, 
precipitation rate is constant with height. Thus, an IWC of 
4gm“ J introduced at 11 km, with the ice morphological and 
size-distribution properties consistent with Airbus measure- 
ments, corresponds to a rain rate of roughly 8 mmh -1 under 
saturated conditions. Figure 12 shows that half as much IWC 
aloft, without a change in ice size distribution or properties, 
corresponds to roughly 4 mmh -1 throughout a saturated col- 
umn. 

Unlike the constant precipitation rate profile in a satu- 
rated column at steady state, the simulated CWC profile de- 
creases monotonically and drastically between 1 1 km and the 
surface. Most dramatically, melting ice causes a rapid in- 
crease in fall speed (consistent with observed MDV), which 
is accompanied by a rapid decrease in CWC to maintain 
a uniform (steady-state) precipitation rate. Above the melt- 
ing level, IWC undergoes a more gentle decrease during sed- 
imentation; simulations that roughly reproduce the trend in 
Z e with height above the melting level also exhibit a roughly 
30-50% decrease in IWC between 11 and 5 km. Thus, ag- 
gregation consistent with increasing Z e corresponds to a re- 
duction of IWC as ice sediments faster. 


6 Meteorological context 

The Australian Bureau of Meteorology C band scanning po- 
larimetric radar at Darwin (C-POL; Keenan et al., 1998) 
offers a wider view of rain rates, which are retrieved at 
2.5 km in altitude following Bringi et al. (2001), and a three- 
dimensional view of hydrometeor class aloft (Keenan, 2003; 
May and Keenan, 2005). Of the two quiescent stratiform 
periods identified in Fig. 8, the second turns out to be far 
more extensive horizontally, with higher peak rain rates seen 
by disdrometer (Fig. 9). Figure 13 shows 2.5 km cross sec- 




Figure 12. Same as Fig. 1 1 except that the initial ice size distribu- 
tion is scaled by a factor of 0.5. 

tions of Z e and rain rate from C-POL measurements near the 
time of maximum stratiform rain rate at the S-band profiler, 
whose location is indicated with an asterisk. Areas outlined 
in black are those identified as convective using the Steiner 
et al. (1995) textural algorithm as in Fridlind et al. (2012a); 
surrounding areas are classified as stratiform. At the time 
shown, the S-band is embedded within an extensive strati- 
form area where most rain rates are 2^4 mmh -1 , although 
higher rain rates of 4-8 mmh -1 can be found locally roughly 
50 km northwest of the profiler. Thus, maximum stratiform 
rain rates in these observations appear roughly consistent 
with the maximum measured IWC aloft, sedimenting under 
saturated conditions, in our calculations, but such rain rates 
are not widespread. The locally stronger stratiform rain rates 
could also be attributable to enhanced precipitation from 
mesoscale ascent concentrated at lower elevations or other 
processes that could enhance rain rate locally by other means, 
rather than high IWC far aloft. 

Although convective rain areas that contain some amount 
of hail are found north and south of the stratiform rain area 
over the S-band in Fig. 13, Z e does not exceed 30dBZ at 
1 1 km anywhere within the C-POL domain at that time. Fig- 
ure 14 shows the C-POL fields at the time of the strongest 
convective rain over the S-band, by contrast. Here the S- 
band is now embedded within a linear band of convection 
containing a continuous elongated region with some amount 
of hail according to C-POL definitions. However, at 1 1 km 
there remain only several tiny areas of 30dBZ Z e ; these ar- 
eas are periodically slightly larger during some periods ear- 
lier and later, but are never extensive during this event (not 
shown). Vertical wind speeds retrieved at 9 km from the pro- 
filers reach 4-10 ms -1 during roughly an hour-long period 
containing convective cells over the S-band, including the 
time of Fig. 14. These would presumably be associated with 
moderate turbulence. 
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Figure 13. Cross sections from C-POL radar fields when the S- 
band experienced stratiform rain at 23.8958 (cf. Fig. 8): Z e and rain 
rate at 2.5 km altitude, Z e at 1 1 km altitude, and column maximum 
hydrometeor class (see legend). The location of the S-band profiler 
is marked with an asterisk. Areas outlined in black are identified 
as convective (see text) except in the 1 1 km Z e field, where they 
enclose areas greater than 30 dBZ (only present in Fig. 14). 

Separation of rainy areas into convective and stratiform ar- 
eas has an origin in the first-order division of mid-latitude 
squall lines into ice particle source regions, where vertical 
winds commonly exceed ice crystal fall speeds (thus trans- 
porting particles to detrainment regions), and ice particle sink 
regions dominated by ice particle sedimentation (Biggerstaff 
and Houze, 1991). Subsequent analyses usefully extended 
this framework to consideration of tropical storms, the dif- 
fering properties of convective and stratiform region energy 
and water budgets, and their unique consequences for atmo- 
spheric circulation (e.g., Houze, 2004; Schumacher et al., 
2004; Zeng et al., 2013). However, the structure of convec- 
tive regions within a tropical MCS present a large variety 
(Mapes and Houze, 1992, cf. their Fig. 15). In the system 
studied here, convection is sometimes organized into linear 
features, but individual convective turrets are also commonly 
scattered within regions identified as stratiform, consistent 
with otherwise quiescent vertical winds and low rain rates. 

7 Conclusions 

Roughly 100 well-documented jet engine power loss and 
damage events described in the literature to date have been 
characterized by high IWC-low-Z e conditions encountered 
almost exclusively in the vicinity of deep convection under 
conditions of weak to moderate turbulence. During 2010- 
2012 Airbus flight tests from Cayenne, Darwin, and Santiago 
carried a new SEA Robust hot-wire probe and an imaging 
nephelometer used to report ice particle size distributions. 
A chief objective of the flight tests was to find and remain 
within the highest IWC regions, within the limits of safety, in 
a manner similar to that described by McNaughton (1959). 


Figure 14. Same as Fig. 13 except when the S-band experienced 
convective rain at 23.7639. 


Whereas that earlier study measured CWC at 0 to —26 °C, 
the Airbus flights primarily sampled at colder temperatures, 
where many such events have occurred (cf. Grzych and Ma- 
son, 2010). 

Similar to CWC measurements reported by McNaughton 
(1959), IWC measurements of roughly 2 gm -3 were not rare 
at any of the three global regions sampled, and maximum 
IWC exceeded twice that value at 10 km length scales (see 
also Grandin et al., 2014). These IWC measurements are 
roughly a factor of 2 greater than measurements from sim- 
ilar cloud regions documented in the scientific literature to 
date. However, as noted above, a large database of such mea- 
surements is not yet available. In addition, the factor-of-2 un- 
certainty generally assigned to IWC measurements by any 
means, in the absence of past test facilities to provide ground 
truth (Strapp et al., 2008; Baumgardner et al., 2011), could 
also account for some of the difference across various direct 
and indirect measurement approaches. 

An examination of the Airbus nephelometer data in the 
regions of the highest IWC at elevations circa 10 km indi- 
cate that mass size distributions were dominated by particles 
in the 100-500 pm range of area-equivalent diameter consis- 
tently across all locations. The mass size distributions also 
consistently exhibited an unexpected self-similarity, or sta- 
bility of shape, over large changes in IWC. Calculations of 
Z e size distributions from the nephelometer data as shown in 
Fig. 4 make apparent that the largest particles present were 
not measured by the nephelometer (limited to D max nomi- 
nally smaller than 2-3 mm), perhaps unrealistically enhanc- 
ing the stability of mass median diameter over the size range 
seen. On the other hand, similar features have also been re- 
ported in previous data sets and analyses. 

Owing to the great uncertainty associated with both IWC 
and size -distribution measurements, here we describe a com- 
parison of the measurements with ground-based profiling 
radar measurements that offer two independent measure- 
ments sensitive to ice morphological and size-distribution 
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properties over the full atmospheric column observed: Z e and 
MDV. We focus on regions of heavy stratiform rain within an 
MCS observed with remote-sensing instruments over Darwin 
on 23 January 2006, where large-scale vertical motions are 
neglected to first order compared with reflectivity-weighted 
ice fall speeds. During heavy stratiform rain periods, Z e is 
less than 30dBZ at typical Airbus sampling elevations of 
11 km and MDV is found to be roughly 1 ms -1 , consistent 
with that calculated from Airbus flight test data. MDV ap- 
pears relatively weakly correlated with Z e at 1 1 km, which is 
consistent with self-similar ice size distributions. However, 
gravity waves with periods of perhaps one -half of the strati- 
form rain duration appear present in vertical wind speed re- 
trievals available up to 9 km, complicating analysis of MDV 
relationship to Z e . 

Column model simulations initialized with ice morpholog- 
ical and particle size distributions consistent with Airbus Ro- 
bust probe and nephelometer measurements are intended to 
test whether the properties of a steady-state stratiform rain 
column under conditions of heavy precipitation appears con- 
sistent with the measurements. If the column is assumed sat- 
urated with respect to ice and water above and below the 
melting level, then the precipitation rate is independent of 
height, and an IWC of 4gm -3 is found to correspond to 
a rain rate of roughly 8mmh -1 . This is roughly twice the 
maximum rain rate observed by disdrometer or retrieved at 
2.5 km in the 23 January event, but polarimetric retrievals 
of rain rate around the profiler did report such rates locally 
within the stratiform region. It is unknown whether such lo- 
cally enhanced rain rates can be attributed to ice aloft or some 
mesoscale ascent, which seems equally possible. Model sim- 
ulations also indicate that in a column that is 10% under- 
saturated in the vicinity of the melting level, precipitation 
can be rapidly reduced by evaporation. Finally, steady-state 
simulations with aggregation sufficient to explain increas- 
ing Z e with decreasing elevation owing to aggregation alone 
also exhibit roughly 30-50% reduced IWC by the melting 
level owing to increasing mass-weighted fall speed. Thus, 
in a steady-state column under such conditions, IWC is ex- 
pected to be higher aloft than at the melting level. In sum- 
mary, the highest Airbus IWC measurements near Darwin 
do appear grossly consistent with the highest stratiform rain 
rates and ice properties observed during the 23 January 2006 
event over Darwin. 

The objective of Part 1 of this two-part study has been 
to vet the Airbus size-distribution measurements to the de- 
gree possible in comparison with past measurements and 
a uniquely relevant remote-sensing data set. A simplified 
modeling framework served to extend the comparison by re- 
lating IWC aloft with rain rates below. Given a lack of gross 
inconsistencies, here we conclude that size distribution and 
IWC measurements appear to give a reasonable representa- 
tion of high IWC-low-Z e ice properties, albeit with large 
uncertainties. The objective of Part 2 is to examine micro- 


physical pathways by which such properties could likely be 
achieved under relevant updraft conditions. 
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